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Abstract 

We solve exactly a generalized Hubbard ring with twisted boundary condi- 
tions. The magnitude of the nearest-neighbor hopping depends on the oc- 
cupations of the sites involved and the term which modifies the number of 
doubly occupied sites t^B = 0. Although ?7-pairing states with off-diagonal 
long-range order are part of the degenerate ground state, the behavior of the 
energy as a function of the twist rules out superconductivity in this limit. 
A small tAB breaks the degeneracy and for moderate repulsive U introduce 
superconducting correlations which lead to "anomalous" flux quantization. 
PACS numbers: 74.20Mn, 71.27.+a, 71.30+h, 71.28.+d 
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One of the most interesting problems of the physics of highly correlated electronic systems 
is the characterization of the metallic, insulating and superconducting phases, as well as the 
transitions among them. Kohn has shown that the Drude weight is the adequate quantity 
to identify the metal-insulator transition (MIT) while Yang introduced the concept of 
off-diagonal long-range order (ODLRO) to characterize the superconducting nature of a 
metallic phase [§. ODLRO in a// relevant low-energy eigenstates implies a periodicity of 
h/2e in the free energy as a function of a magnetic flux threading a system with annular 
topology, which is referred to as "anomalous" flux quantization (AFQ) In other words, 
AFQ in the ground-state (GS) means -|- vr) = where E is the GS energy and $ 

is the twist angle. It is a necessary but not sufficient condition for superconductivity 0. 
Since ~ d^E/d^^, the function £'($) gives crucial information about the metallic and 
superconducting character of the system |jl],§J^. 

Exactly solvable highly correlated models displaying a MIT or ODLRO are good labora- 
tories to investigate the nature of the MIT and electronic mechanisms of superconductivity. 
The Bethe ansatz solution with twisted boundary conditions (i.e. arbitrary flux) of the one 
dimensional (ID) Hubbard model [0, allowed to apply Kohn's ideas to the MIT in this 
model Very few exact results exist, for electronic models exhibiting superconductiv- 

ity (or dominant superconducting correlations at long distances in ID). Several of them 
are related with the so called //-pairing mechanism, which allows to construct eigenstates 
with ODLRO [p|-[T^. In particular, the widely studied [|1^ Hubbard model with correlated 



hopping. 



H = Hu + Ht = U^ni^^riii + ^ (cJaCja + h.c){tAA (1 - ni„){l - rij^) + 



<tJ>(T 



tBB rii^Uj^ + tAB K<t(1 - nj^) + nj^{l - Ui^)]}, (1) 



has been exactly solved recently in ID in the limit tAB = for open |[10| and periodic p!l|JT2 
boundary conditions. It has been shown that, due to an SU(2) pseudospin symmetry, re- 
pairing states with ODLRO are part of the degenerate GS for moderate on-site repulsion 
U and arbitrary band filling. Unfortunately, the function E{^) has not been obtained and 



then, the AFQ and Dc were not studied. Our main interest in this study is motivated by 
the following two facts: first, the superconducting character of the degenerate GS is not 
obvious, even when states with ODLRO are part of the GS manifold. Second, the GS was 
found to be a Mott insulator for U > Umi = 2D(|t^yi| + l^ssl) at half-filling (density of 
particles n = 1), in a simple cubic lattice in D dimensions, with a MIT for D> 1 |T^,|TB[. 
Strictly in D= 1, however, we find that for n = 1, Dc = \fU, in spite of a vanishing charge 
gap for U < Umi- The possibility of an insulating phase with this feature was first remarked 
by Kohn and we think that this is, to our knowledge, the first non trivial realization of that 
kind of insulators. 

In this Letter, we solve exactly the model (HD with t^B = 0, for twisted boundary 
conditions $1 ($1) for spin up (down) fermions. This allows us to calculate the Drude 
weight Dc and the spin stiffness D^ and to discuss the nature of the MIT as n — 1. The 
behavior of i?($|,<l>j), rules out superconductivity in the model for Iab = 0, at least in 
the one-dimensional case we study here. In addition, we go one step further and show how 
the GS degeneracy is broken in favor of a state with dominant superconducting correlations 
when a finite t^B is allowed, for moderate repulsive U. 

We first consider the Hamiltonian (0) with — t^A = tBB = t > 0, t^B = 0. The other 
possible choices of the sign of tAA and tsB lead to an equivalent model At each site 
i, we introduce two fermions /jo- and two bosons bi„/, where 6^+ = e (empty) and 6j_ = d 
(doublon). The fermions (bosons) transform according to an irreducible representation of 
the local spin (pseudospin) local SU(2) symmetry that Ht possesses for open boundary 
conditions [0. In this representation, Ht of ([l|) in a ring of L sites with twists $0- for 
particles with spin a reads 

Ht = J2H,,+, + Hl,i = -t J2 ifUiJiabl'h+ia' + h.c.) 

1=1 i=l,acr' 

-tT^lflfLM^-blbi^^ + e-*--6l_6L-) + h.c.]. (2) 

The numbers = J2ifLfi<T, Ne = J2i4a(^i<T, Nd = T.id\^di„ are all conserved 
when tAB = 0. In each subspace with fixed N-^, N^, N^, N^, any state has the form 
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rimLi n^i //(j)<TO)|0)' where j labels the Nf = + Ni fermions from left to right 

and i{j),a{j) denote the position and the spin of the jth fermion. Similarly 1(171) , a' {m) 



(with i(m + 1) > i{m)) are the position and the pseudospin |10| of the mth boson. The 
number of bosons is Ni, = N^. + = L — Nf. Because of the completeness relation 
J2a ficrfi^T + blf^bia = 1, {^(j)} and {i{m)} are complementary sets. 

For the periodic case (^^ = 0), the Hamiltonian is invariant under cyclic permutations 
of the fermions and bosons and it is convenient to work in the basis of the irreducible 
representation of the direct product group Cat^ (8) Cat,, (T^- ^^^^ appropriate 
weighted representations |^ to cancel out the difference in phases in Hi l in order to map 



the problem into one of spinless fermions with twisted boundary conditions. With this 
objective in mind, we think the ring as a periodic system in which //^^i = fjp ej_,_^ = e|, 
but fi^i^i = e~*^*T-*i) J^t^^ c?|_|_^ = e~*(*T+*i)c?|. Using these boundary conditions, it can be 
verified that 

Hl,i = e'^^Y.fUjLA-.L.'h.', (3) 

acr' 

as we want. 

We look for a basis of many particle states transforming as irreducible representations of 
^Nf ^^Nt under the above mentioned boundary conditions. The part of these states which 
describes the singly occupied sites can be constructed using the operators 

F\m},{ki}) = uhrI[fL 

ft -^Ve-^^ft , _2vn^,-$,) 
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= ^ E e"^'V4-)T(i - fhnkM)^ (4) 



where in contrast to the wave numbers ki, the ki are not shifted [ki Nf/2TT is integer), and 
the set of ki is fixed and chosen in such a way that J2i ki= (tt) for Nf odd (even). The ui are 
N^ different integers lying in the interval [0,Nf — 1], and each of the Nfl/lN-^lN^l) possible 



choices of the set of ui define a spin configuration. It is easy to see that under cychc permuta- 
tion Ctv^, which carries each fermionic position to the right C^^.F'^ = —{—l)'^f exp(i J2i ki)F^ . 
In a similar way, using a transformation that interchanges spin and pseudospin ||21|, the 
pseudospin configuration can be described by an operator B'^ {{i{m)}, {k'l}), such that 
Cjv,5t = exp{iY.ik'i)B^, with the A^'^ different k'l = [2nu' - ($1 + $|)]/iVb. (details 
will be given elsewhere). The (non-orthonormal) basis states that we use are denoted by 
\mj)}Ak},{k'})=B^F^O). 

Ht permutes a fermion and a nearest-neighbor boson. The cyclic orders of fermions 
and bosons are conserved. Thus, the numbers {k} and {k'} are conserved. We drop these 
indices for simplicity. Hii^i\tp{i{j)}) = unless one and only one of the sites / and / + 1 
is contained in We restrict ourselves to these states in the following discussion. It 

can be seen that for / < L, = — where {«'(/)} differs form 

in the position of one fermion only, which is shifted from site / to / + 1 or conversely. 
If tiNj) = L, then ifi,i|vl/({^(j)})) = -te^*TC^^Q;|^{^'(j)}) = t(-l)^/ exp[z(E£ fc, - 
k'l + ^^)]\ip{i'{j)}), where = 1, and for j < Nf, + 1) = When A^^ = A^^ = 
0, Ht takes the form of a problem of spinless fermions with flux The above equations 
show that in the general case, the problem takes the same form, with an effective flux: 



The energy of the system is given by 

E = ~2tY: cosi ^'"'"'"''^'" ) + UNd, (6) 
3=1 ^ 

where the Nf integer numbers ^''j should be different and can be chosen in the interval 
[—L/2 + 1, L/2]. When A*";, = $| = $| = 0, the energy takes a similar form as that derived 



by Ogata and Shiba I^Sj for the infinite U Hubbard model from the Bethe ansatz equations. 
Note, however, that our basis states do not correspond to a definite total spin in general. 
Also, our basis is complete, and describes the energy of any state, while the regular 



Bethe ansatz states do not form a complete basis ||23 



We are in position now to examine the ground state energy as a function of the fluxes. 
For fixed Nf = UfL and total number of particles N = nL = Nj + 2Nh, minimization of 
leads to 



sm(raj7r) 



E,(<l>t,<fi) = UN, - 2t— ^cos(¥.), (7) 

where ip = ^eff/L for Nf odd and (f = {^eff ~'^)/L for Nf even. The value of U determines 
Nd for the GS. For each $|, $j, the numbers N^, Ni, as well as {u}, {u'} should be chosen to 
minimize \(p\ (module 2tt). It can be easily seen that in the simplest case $1 = $1 = and 
U = that the GS is highly degenerate (many choices of quantum numbers lead to ip/27r 



integer). For U > Uc = — 4tcos(7™), it was found |1T0|-|T2| that the double occupation is 
forbidden in the GS (A^^ = 0), and we recover the solution of the U = +00 Hubbard model 
with twisted boundary conditions. In this case _E'($| + 27r/N,^^ + 2-k/N) = £'(<l>|,<l>|), 
since the shift in $o- can be absorbed decreasing one of the z/ by 1, what is always possible 
if 7^ A^^l 7^ Nf. For $| = $j, this result has been obtained previously ||2^. Similarly, 
for the more interesting case with A*"^ 7^ 0, a change in both $0- by 271 /\Nh — N^l can be 
counterbalanced by a change in the {i^'}, leading to 

^.(^T + + ^^3^) = E,{<l>,, (8) 

for L N, whereas for a half-filled system Eg depends only on the difference $| — a 
behavior typical of an insulator. For U < — 4cos(7r?2), A^^^ > and r^-pairing states with 
ODLRO are present in the GS ||10|-[T^. However, we do not find AFQ, but a periodicity 
which depends on the particle-density n. An example for finite chains is shown in Figs. 1 
(a-b). The number of peaks of -Eg($t = "J^j, = *^') for < $ < 27r is at least L\n — 1|, 
diverging in the thermodynamic limit, while the height of each peak decreases as as 
in the U = +00 Hubbard model. The response of the system to the flux is like that of a 
single particle with charge L — N or larger. One might ask whether a collection of weakly 
coupled chains behaves like a superfiuid of these particles. However, since the compressibility 
diverges in the interesting regime |]lOl , charge can be transferred between chains without cost 
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of energy, and the response to the flux of different chains does not add coherently. Thus, 
the system does not show the Meissner effect [0-^. We should also note that the SU(2) 
?7-symmetry which allows for the construction of eigenstates with ODLRO is broken in the 
presence of a flux. It can be seen from that a twist of the form = — couples with 
the spin degrees of freedom in the same fashion as a twist $| = $| affects the pseudospin 
ones. While in the first case, the total spin invariance is broken, remaining only as a 
good quantum number, the total pseudospin invariance is broken in the second case, with 
7]^ = 1/2(L — A^) fixed by U and the chemical potential. Thus, the ?7-pairing states do not 
necessarily give rise to superconducting currents in the presence of an external flux. The 
physics in the region with U < ~4:t, where nf = [lTO|-p!^ is more obvious. In this case. 



there are also r^-paired states with ODLRO in the degenerate GS. However, these states are 
static, and from (|^) = UN^ for all $0-. This result might been anticipated from 

the form of Ht. This demonstrates that the ODLRO of the ?7-paired states is not a sufficient 
condition for the existence of superconductivity. This fact has not been noted in previous 
related works. 

The computations of Drude weight and the spin stiffness lead to 

^e=-^^|.^o = -(^3^) sm(vrn,) 
^ L -$) , t,n-r-n\.^ . , , 



where is a function of U /t and n that can be obtained minimizing (|^) [|T0|,[T2|. For half 



filling, D(. = 0, Vt/. The result is not surprising for \U\ > 4t, where the system is a Mott 
insulator, but rather unexpected otherwise. The symmetry between spin and pseudospin 
degrees of freedom [^], becomes explicit by replacing Ue — Ud = 1 — n, rib = 1 — nf in the 
expression of D^. Dg vanishes in the sector S*^ = and with it, the inverse of the magnetic 
susceptibility 10, as a consequence of the spin degeneracy. Analogously, for \U\ < 4t, not 
only Dc vanishes at half-filling, but also the inverse of the charge compressibility. The 
system is an insulator, in spite of a vanishing charge gap. The behavior of the Drude weight 
for n — i> 1 is the same as that of a system of carriers with effective mass diverging as 



(l-n;)V(l-n2). 

The particular behavior of the energy as a function of flux, is a consequence of the 
existence of excitations associated to the fermionic charge, spin and pseudospin degrees of 
freedom, which is explicit in Eqs. @ and (|^). This is similar to the case of the infinite 
U Hubbard model where spin and charge decouple. The ground state is highly degenerate 
as a consequence of the rich symmetry structure of Eq. (|l]) when — \tBB\ = tAB = 0. 
In particular, for an open chain there is a local spin and pseudospin symmetry at each 
site Q. Since Iab = is an accident rather than a generic feature of any one-band 
model, it is very important to discuss the effect of a finite Iab, particularly taking into 
account that this term hfts the GS degeneracy. When tAB 7^ 0, the sign of tAB or those 



of tAA and tBB simultaneously, can be changed using symmetry properties [|T^,|T5|, but 
models with different tAA^BB are not equivalent. In the following, we consider the case 
tAA = tBB = —t, which interpolates between two exactly solvable cases: the one considered 
above and the Hubbard model. This case preserves the SU(2) pseudospin symmetry when 
tAB 7^ |T^JT^. Previous numerical studies suggest that the ID system for small U and 
n ~ 1 is a Tomonaga-Luttinger liquid with dominant superconducting correlations at large 
distances |T^|T^. We have studied numerically = $| = $) for finite systems with 



fixed densities n ^ 1, Ue^Ud 7^ at tAB 7^ 0. An example is shown in Fig. 1 (c). We 
find that a small tAB gives rise to AFQ in the finite systems. The cusp near $ = 7r/2 
decreases with increasing length of chain. This behavior is typical of a system with power 
law superconducting correlations at long distances rather than a state with true ODLRO. 

What is the origin of the superconducting correlations? Is it related with the r^-pairing? 
To answer these questions let us begin by noting that the (nondegenerate) GS for tAB — ^ 
is exactly known in the case N^Nh = (for U > — 4tcos(7™) and U > p!0|-p]2|). In 



this case, the low energy physics of the model becomes equivalent to that of a Hubbard 
model with interaction Uh = f^U/t^^ 00 (as can be easily seen eliminating in (^ the 
term in tAB through the standard canonical transformation). In this limit Ogata and Shiba 



2^ have shown that the GS wave function can be factorized in two parts: one describing 



the position and the other the spin of the Nf fermions. The first factor corresponds to 
the GS of a Heisenberg chain with Nf sites. This GS wave function can be mapped into 
the corresponding one for f/ < and a magnetic field high enough to ensure N-^Ni = 
using the transformation that interchanges spin and pseudospin It is natural to 

expect that the effect of a small tAB is to introduce antiferromagnetic correlations between 
spins and pseudospins in the general case. A straightforward generalization of the above 
analyzed two cases, led us to propose an ansatz for the GS in the limit Iab consisting of 
three factors, describing the positions of fermions and bosons and the spin and pseudospin 
variables. The first two factors are those of the GS Bethe ansatz solution of the U = +oo 
Hubbard model. The last one is the GS of a Heisenberg model for the pseudospin variables, 
which is also the GS of the large negative U Hubbard model in a system with A*";, sites 
and 2Nii particles ||2^. We have computed the overlap of our ansatz with the exact GS 



obtained from exact diagonalization in different chains (up to L = 12 sites), and we found 
that it is equal to (1 — a^t^^)^/^, with a ~ 2 — 4, for 7^ t^B < 0.1, confirming our 
conjecture for t^B ~^ 0. It is easy to verify with our ansatz that in the thermodynamic 
limit, for L ^ 00, the pair correlation function C{1) = (c|_^;|c|_,_^Cj|Cj|) can be expressed in 
terms of the corresponding correlation function of the large \U\ attractive Hubbard model 
Ch{1) for density (n — nf)/{l — Uf) as Cuil) = (1 — nf)[CH{l')]av, where the average V 
is centered around L(l — nj). Thus, the superconducting properties of the system are 
essentially those of the large U attractive Hubbard model with dilute superfluid density. The 
superconducting properties of the model are not related with the ?7-pairing. For Iaa = ^bb-, 
the generators of the total pseudospin algebra are 77" = Z]i(~l)*4t4i? 'n^ — iv V ^^^^ 
T]^ = (l/2)X]j(l — Z^o-^io-)- The r^-pairing mechanism applies rj' to an eigenstate with 
?7 7^ to obtain eigenstates with rj^ < t] which possess ODLRO p|-p!^j. However, since by 
construction, our ansatz for the GS has 1]^ = rj, and the true (non-degenerate) ground state 
has the same quantum numbers, it cannot be the result of applying 77" to any eigenstate. 
Exact diagonalization results show that this is also the case in 2D, even for tAB = | |17|| . 
In summary, we have shown that at least in the ID generalized Hubbard model (|l]) for 
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I^aaI — \tBB\ = tAB = 0, the r^-pairing does not lead to superconductivity. The possibihty 
of constructing eigenstates with ODLRO using the SU(2) symmetry does not guaranty the 
existence of superconducting currents giving rise anomalous flux quantization and Meissner 
effect. The ODLRO must be analyzed in the presence of a finite magnetic flux threading 
the ring. This fundamental fact is in the spirit of the proposals of Refs. [0-^. We have 
also examined the character of the metal-insulator transition near half fllling and we have 
presented strong evidence that the GS degeneracy is broken in favor of a GS with dominant 
superconducting correlations in ID when a small Iab is turned on. 
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Figure Captions 

1. Ground state energy as a function of twist angle (a) for {IaaI = l^ssl = 1, ^ab = 0, 
density n = 2/3 and U = 0; (b) same as (a) with U > 4; (c) for tAA = tsB = —1, tAB — 
-0.2, C/ = 0, n = 2/3 
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